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Abstract 

We present the first observation of r lepton decays to hadronic final states with 
a (/)-meson. This analysis is based on 401 fb~^ of data accumulated at the Belle 
experiment. The branching fraction obtained is B{t~ — s- cjiK^Vj-) = (4.05 ± 0.25 ± 
0.26) X 10"^ 

Key words: tau, phi 
PACS: 13.35.Dx, 14.40.Cs 



1 Introduction 



Hadronic r decays witli a meson in the final state are valuable to investigate 
QCD at a low mass scale. However, they have never been observed due to their 
small branching fractions. The decay r~ (pK'Vj. is Cabibbo-suppressed and 
further restricted by its small phase space, while the decay r" (pTi'v^. is sup- 
pressed by the OZI rule although it is Cabibbo-allowed (Fig. 1). The branch- 
ing fraction of the former can roughly be estimated by scaling the analogous 
Cabibbo-allowed decay r~ K*K~Ut^ [1] by tan^ 9c and the ratio of the phase 
space of the two decays, resulting in B{t~ — ^ cjyK'VT-) ~ 2x 10~^. Similarly, the 
vector dominance model predicts i3(r~ ^vr"!^^) = (1.20 ± 0.48) x 10^^ [2], 
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whereas the CVC upper hmit following from the cross section for e~^e ^vr' 
is B{t- (pn-Ur) < 3 X 10"^ at the 90% confidence level [3]. 




Fig. 1. Diagrams for r (j)K (left) and r <j)7r (right). 

Previously, the CLEO collaboration searched for these decays using 3.1 fb~^ 
of data taken on the T(4S') resonance. They set upper limits of i3(r~ 
(t)K"Ur) < (5.4 - 6.7) X 10^5 and B{t- (pTx~Ur) < (1.2 - 2.0) x 10^^ at the 
90% confidence level, depending on the mechanism assumed for the decay [4]. 
Here we report the first measurement of the r~ — > (j)K~Ur decay. (Throughout 
this paper charge-conjugate states are implied.) We also observe for the first 
time the decay r~ — > 07r~z/^, but it is treated here as a background process, 
together with the kinematically allowed but phase-space suppressed decays 
T~ 07r^(?T,7r)z/^ ^ n < 4). The result is based on a data sample of 401 
fb~^ corresponding to 3.6 x 10® t~^t~ pairs collected near the T(45') resonance 
with the Belle detector at the KEKB asymmetric-energy e~^e~ (3.5 on 8 GeV) 
collider [5]. 

The Belle detector is a large- solid- angle magnetic spectrometer that consists of 
a silicon vertex detector, a 50-layer central drift chamber, an array of aerogel 
threshold Cerenkov counters, a barrel-like arrangement of time-of-fiight scin- 
tillation counters, and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals located inside a superconducting solenoid coil that provides a 1.5 T 
magnetic field. An iron fiux-return located outside of the coil is instrumented 
to detect Kl mesons and identify muons. The detector is described in detail 
elsewhere [6]. Two inner detector configurations were used. A 2.0 cm radius 
beampipe and a 3-layer silicon vertex detector were used for the first sample 
of 158 fb~^, while a 1.5 cm radius beampipe, a 4-layer silicon detector and a 
small-cell inner drift chamber were used to record the remaining 244 fb~^ [7]. 



2 Event Selection 



We look for r (pK candidates in the reaction e~^e t^t with the 
following signature: 

^signal ^(p + K^ + (missing) 
^ K+K~ 
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^ttg We)"^ + n{< 1)7 + (missing), 

where 'missing' denotes other possible daughters not reconstructed. The de- 
tection of mesons rehes on the — >• K^K~ decay {E = (49.2 ± 0.6)% [1]); 
the final evaluation of the signal yield is carried out using the K~^K^ invariant 
mass distribution. 

The selection criteria described below are determined from studies of Monte- 
Carlo (MC) simulated events. The background samples consist of t~^t~ (1570 
fb~^, which does not include any decay mode with a meson) and qq contin- 
uum, B^B^, B^B^ and two-photon processes. For signal, we generate samples 
with 2 X 10^ (pR'^T, (pR'TT^UT, (f)TX~Vj. and 07r~7r°i^^ events. 

The transverse momentum for a charged track is required to be larger than 
0.06 GeV/c in the barrel region (—0.6235 < cos^ < 0.8332, where 9 is the polar 
angle relative to the direction opposite to that of the incident e"*" beam in the 
laboratory frame) and 0.1 GeV/c in the endcap region (—0.8660 < cos^ < 
—0.6235, and 0.8332 < cos 6' < 0.9563). The energy of photon candidates is 
required to be larger than 0.1 GeV in both regions. 

To select a r-pair sample, we require four charged tracks in an event with 
zero net charge, and a total energy of charged tracks and photons in the 
center-of-mass (CM) frame less than 11 GeV. We also require that the missing 
momentum in the laboratory frame be greater than 0.1 GeV/c, and that its 
direction be within the detector acceptance, where the missing momentum is 
defined as the difference between the momentum of the initial e^e~ system, 
and the sum of the observed momentum vectors. The event is subdivided 
into 3-prong and 1-prong hemispheres according to the thrust axis in the CM 
frame. These are referred to as the signal and tag side, respectively. We allow 
at most one photon on the tag side to account for initial state radiation, while 
requiring no extra photons on the signal side to reduce the qq backgrounds. 

We require cos < —0.6 to reduce backgrounds from other r decays 

and qq processes, where ^'^jij^g^.j^jigg is the opening angle between the thrust axis 
(on the signal side) and the missing momentum in the CM frame. In order to 
remove the qq background, we require that the invariant mass of the particles 
on the tag side (if a 7 is present) be less than 1.8 GeV/c^ (^ 771^). Similarly, 
the effective mass of the signal side must be less than 1.8 GeV/c^. Moreover, 
we require that the lepton likelihood ratio P^/e be greater than 0.1 for the 
charged track on the tag side. Here P.j. is the likelihood ratio for a charged 
particle of type x {x = ^, e, K or tt), defined as = L^/ {J2x ^x), where is 
the likelihood for particle type hypothesis x, determined from responses of the 
relevant detectors [8]. The efficiencies for muon and electron identification are 
92% for momenta larger than 1.0 GeV/c and 94% for momenta larger than 
0.5 GeV/c, respectively. 
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We require that both kaon daughters of the candidate have kaon hkehhood 
ratios Pk > 0.8 and cos^ > —0.6. The kaon identification efficiency is 82%. 
To suppress combinatorial backgrounds from other r decays and qq processes, 
we require that the momentum be greater than 1.5 GeV/c in the CM frame. 
After these requirements, the remaining contributions from B^B^, B'^B', 
Bhabha, fi pair and two-photon backgrounds are neghgible. 

To separate (pK'Vr from (pTT'Ur, the remaining charged track is required to 
satisfy the same kaon identification criteria as the (p daughters. The r+r^ 
and qq backgrounds are reduced by requiring that the opening angle (O^k) 
between the and K~ in the CM frame satisfy cos^^^ > 0.92, and that the 
CM momentum of the (j)K~ system be greater than 3.5 GeV/c. For (j)'ir~i>r, we 
require that the charged track be identified as a pion, > 0.8, and that the 
opening angle between the and n~ in the CM frame satisfy cos^^^ < 0.98. 
This last requirement suppresses the background from (pK^Ur and 

Figure 2(a) shows the K^K~ invariant mass distribution after all r~ 
(f)K~Ur selection requirements. As there are two possible K^K~ combina- 
tions from the K~K^K~ tracks on the signal side, this distribution has two 
entries per event. Therefore, the signal MC shape includes a long tail due to 
the wrong K^K~ combination. Non-resonant backgrounds arise mainly from 
T~ K~t\:~Vt-, which has a branching fraction of i3 = (1.53 ± 0.10) x 

10^'^ [1]. Small contributions are expected from qq processes as described be- 
low. 



3 Signal and background evaluation 

The detection efficiencies e for r~ (pK'Ur and the cross- feed rates from 
(pK'TT^Vr and 07r~z/T- are evaluated, as listed in Table 1, from MC simulation 
using KKMC [9], where the V — A interaction is assumed at the vertices and 
the final hadrons decay according to non-resonant phase space. The efficiencies 
include the branching fraction for K^K~ . 

Table 1 

Detection efficiencies e and cross-feed rates (%), from MC simulation. The errors 
are from the MC statistics. 







Decay modes 




Candidates 




(pTTl/ 




r (j)Kv 


1.826±0.009 


0.049± 0.002 


0.328±0.006 




0.110±0.002 


1.663± 0.014 


0.009±0.001 



The signal yields are extracted by a fit to the K^K invariant mass distribu- 
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M(K,K) GeV/c^ 

Fig. 2. K^K~ invariant mass distributions for (a) t~ (pK'v-T and (b) 

(pTT^Ur- Points with error bars indicate the data. The shaded histograms 
show the expectations from t~^t~ and qq background MC simulations. The open 
histogram is the signal MC with B{t~ —>■ (f)K~v.r) = 4 x 10~^ in (a) and 
B{t~ — >• (l)'K~Vr) = 6 X 10~^ in (b). The curves show the best fit results, and 
the dashed curves indicate the non-resonant background contributions. See the text 
for details. 

tion. For signal, we use a p-wave Breit-Wigner (BW) distribution convolved 
with a Gaussian function (of width a) to account for the detector resolution. 
The (j) width is fixed to be = 4.26 MeV/c^ [1] but a is allowed to float. First- 
and second-order polynomial background functions are used for r~ (pK'Ur 
and (pTT'i^T- decays, respectively. The fit results are also shown in Fig. 2. The 
obtained signal yields are N^ku = 573 ±32 and N^nu = 753 ±84. The a's from 
the fits are 1.2 ± 0.3 MeV/c^ and 1.2 ± 0.7 MeV/c^ for (pR-u^ and (pir-u^, 
respectively, which are consistent with MC simulation. 

MC studies show that only the (pTT"!/^-, r~ (f)K~TX^Vj. and qq samples 

yield significant contributions peaking at the mass. The contributions of 
other backgrounds are less than 0.01% and can be neglected. The contribution 
of T~ </)7r~z/^ events to the (pK'Vj. sample is estimated using N^^^y and the 
misidentification rate, as discussed below. Other contributions are estimated 
as follows. 

To evaluate the branching fraction and background contribution from r~ 
(j)K~7r'^h'T-, we select vr^ 77 candidates and combine them with (pK^Vj. 
combinations that satisfy the requirements listed above. The signal yield is 
estimated by fitting the resulting K^K~ invariant mass distribution with a 
p-wave BW distribution plus a linear background function, as shown in Fig. 3. 
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The resulting yield is 8.2 ± 3.8 (pKir^v events. Using a detection efficiency 
^(I>k-kOv = (0.396±0.007)% obtained from MC simulation, and an e+e~ t^t~ 
sample normalization A^,-^ = 401 fb~^ x 0.892 nb = 3.58 x 10*, we obtain a 
branching fraction B{t~ -h- (f)K~TT^UT-) = (2.9 ± 1.3) x 10^^. However, this 
must be corrected for the unknown contamination of r~ (j)'K~'K^{n'K)vr (0 < 
n < 3) decays. Using this value, we estimate the r" (pK^n^Ur background 
in the r" (pK'Ur sample to be N^k-k^v = (6.8 ± 3.1) events, given a cross- 
feed rate for r" (pR-n^u^ to the r" ^K'u^ sample of (0.328 ± 0.006)% 
(see Table 1). 




1.03 1.04 1.05 
M(K,K) GeV/c^ 

Fig. 3. K^K~ invariant mass distributions for r~ — > cpK' Points with error 

bars indicate the data. Histograms show the MC expectations of r-pairs (shaded) 
and signal (open) with a branching fraction of 3 x 10~^. The sohd curve shows the 
best fit result and the dashed curve shows the non-resonant background contribu- 
tion. 



From a MC study, we find a qq contamination of Nqq = 6.6 ± 2.5. To take 
into account the uncertainty in production in the qq MC, we compare MC 
results with enriched qq data by demanding that the effective mass of the 
tag side be larger than 1.8 GeV/c^. With this selection, the background is 
qq dominated and the other backgrounds are negligible. The yield in data is 
262±21 events, and the yield in the qq MC is 117±10 events. We subsequently 
scale the above qq background estimate by the ratio / = 2.23±0.26; the result 
is Naa = 14.8 ± 5.8 events. 



4 Results 



The peaking backgrounds described above, r~ — > (pK'n^Vr and qq, are sub- 
tracted from the signal yield, leaving A^'^xi^ = (573±32) — (6.8±3.1) — (14. 8± 
5.8) = 551 ±33 events. To take into account cross-feed between r~ (pK'Ur 
and T~ (pir'Ur due to particle misidentification {K ^ vr), we solve the 
following simultaneous equations: 

N^Ku = 2Nrr [e^Ku X B^Ku + 4nu >< ^^^u) , (1) 
NffjTri^ = 2Nt-t (e^Xy ^ ^^Ku + e(/.7r!/ X B<PTTuj , (2) 
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where B^ku and Bfp^i, are the branching fractions for r~ ipK'Vr and r~ 
07r~z/T-, respectively. The detection efficiencies, e's, are hsted in Table 1. The 
factor e^^j^ is the efficiency for reconstructing r" — > (pK'Uj. as r~ ^ 07r~i/^ 

while e^^J' is the efficiency for reconstructing r~ ^ (pix^VT as ^ (pK^u-r. 
The resulting branching fraction for r~ — (pK^u^ is 

^^K. = (4.05 ± 0.25) X 10-^ (3) 

where the uncertainty is due to the statistical uncertainty in the N^ku and 
Nif)T^^ terms. The uncertainty in the detection efficiencies, e's, will be taken into 
account in the systematic error. The result for Btpnu is B^t^u = (6.05 ± 0.71) x 
10~^; however, small background from r~ —>■ (j)7r~{n7r)i>r {I < n < 4) decays 
is included and must be subtracted to obtain the final branching fraction. 

The systematic uncertainties are estimated as follows: The uncertainties in 
the integrated luminosity, t~^t~ cross-section and trigger efficiency are 1.4%, 
1.3% and 1.1%, respectively. Track finding efficiency has an uncertainty of 
4.0%. Uncertainties in lepton and kaon identification efficiencies and fake rate 
are evaluated, respectively, to be 3.2% and 3.1% by averaging the estimated 
uncertainties depending on momentum and polar angle of each charged track. 
To evaluate the systematic uncertainty of fixing in the BW fit, we calculate 
the change in the signal yield when is varied by ±0.05 MeV/c^ (the uncer- 
tainty quoted by the PDG) [1]: the result is 0.2%. The branching fraction for 
(j) — i> K^K~ gives an uncertainty of 1.2% [1]. The signal detection efficiency 
e^pKu has an uncertainty of 0.5% due to MC statistics. A total systematic un- 
certainty of 6.5% is obtained by adding all uncertainties in quadrature. The 
resulting branching fraction is then 

B{t~ (pK-Ur) = (4.05 ± 0.25 ± 0.26) x 10"^ (4) 

Finally, we consider the possibility that a resonant state contributes to the 
final (f)K~ hadronic system. We generate a resonant MC with the KKMC sim- 
ulation program. The weak current is generated with a.V — A form while the 
(f)K~ system is assumed to be produced from a 2-body decay of a resonance. 
In Fig. 4(a), the (j)K~ mass distribution for data is compared to MC; the 
combinatorial background is subtracted using the K^K~ sideband. The MC 
distributions correspond to (M, F) =(1650, 100) MeV/c^, (M, F) =(1570, 150) 
MeV/c^, and also non-resonant phase space. Figure 4(b) shows the </)'s angu- 
lar distribution in the (pK' rest frame (cos a) , where the negative of the lab 
frame direction in the (pK' frame is taken as the reference axis. It indicates 
an isotropic distribution in the (pK' system. For both the invariant mass and 
angular distributions of the (j)K~ system, the phase space MC reproduces the 
signal distribution well. We therefore neglect systematic uncertainty due to 
possible resonant structure. On the other hand, the 1650 MeV/c^ state as- 
sumed in the CLEO search [4], indicated by the dotted histogram in Fig. 4(a), 
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clearly cannot account for the entire signal. If production via a single reso- 
nance is assumed, the best agreement with data is found for a mass and a 
width of ~1570 MeV/c^ and ~150 MeV/c^, respectively, as shown by the dot- 
dashed histogram. However, since the shape of the resonant MC is similar to 
the phase-space-distributed MC, we cannot draw any strong conclusions about 
an intermediate resonance with F ~ O(100MeV/c^) in this narrow mass range 
of -250 MeV/c2. 

I 50 
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::4o 

in 

I 30 

(0 

I 20 
E 

10 
O 
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Fig. 4. (a) invariant mass and (b) angular distributions for the (l)K~ system. The 
non-(/)-resonant backgrounds are subtracted using the sideband spectra. Points with 
error bars indicate the data. The open histogram shows the phase space distributed 
signal MC, and dotted and dot-dashed histograms indicate the signal MC mediated 
by a resonance with M = 1650 MeV/c^ and F = 100 MeV/c^ and M = 1570 
MeV/c^ and F = 150 MeV/c^, respectively. In the MC, a branching fraction of 

4 X 10~^ is assumed, (b) (/>'s angular distribution in the (l)K~ rest frame, where the 
negative of the lab frame direction in the (pK^ frame is taken as the reference axis. 

5 Conclusion 

Using 401 fb~^ of data, we make the first observation of the rare decay r~ 
(f)K~Vr- The measured branching fraction is 

B{t^ (pK^Ur) = (4.05 ± 0.25 ± 0.26) x 10^1 (5) 



Acknowledgments 





10 



We gratefully acknowledge the essential contributions of Mari Kitayabu, which 
are described in her bachelor thesis at Nagoya University. We thank the KEKB 
group for the excellent operation of the accelerator, the KEK cryogenics group 
for the efficient operation of the solenoid, and the KEK computer group and 
the National Institute of Informatics for valuable computing and Super-SINET 
network support. We acknowledge support from the Ministry of Education, 
Culture, Sports, Science, and Technology of Japan and the Japan Society for 
the Promotion of Science; the Austrahan Research Council and the Australian 
Department of Education, Science and Training; the National Science Founda- 
tion of China and the Knowledge Innovation Program of the Chinese Academy 
of Sciences under contract No. 10575109 and lHEP-U-503; the Department of 
Science and Technology of India; the BK21 program of the Ministry of Edu- 
cation of Korea, the CHEP SRC program and Basic Research program (grant 
No. ROl-2005-000-10089-0) of the Korea Science and Engineering Founda- 
tion, and the Pure Basic Research Group program of the Korea Research 
Foundation; the Polish State Committee for Scientific Research; the Ministry 
of Science and Technology of the Russian Federation; the Slovenian Research 
Agency; the Swiss National Science Foundation; the National Science Council 
and the Ministry of Education of Taiwan; and the U.S. Department of Energy. 

References 

[1] W.-M. Yao et al. (Particle Data Group), J. Phys. G 33, 1 (2006). 

[2] G. Lopez Castro and D. A. Lopez Falcon, Phys. Rev. D 54, 4400 (1996). 

[3] S. L Eidelman and V. N. Ivanchenko, Nucl. Phys. B (Proc. Suppl.) 55C, 181 
(1997). 

[4] P. Avery et al. (CLEO Collaboration), Phys. Rev. D 55, R1119 (1997). 

[5] S. Kurokawa and E. Kikutani, Nucl. lustrum. Meth. A 499, 1 (2003), and other 
papers included in this Volume. 

[6] A. Abashian et al. (Belle Collaboration), Nucl. Instr. and Meth. A 479, 117 
(2002). 

[7] Y. Ushiroda (Belle SVD2 Group), Nucl. Instr. and Meth. A 511, 6 (2003). 

[8] K. Hanagaki et al, Nucl. Instr. and Meth. A 485, 490 (2002); A. Abashian et 
a/., Nucl. Instr. and Meth. A 491, 69 (2002). 

[9] S. Jadach, B.F.L. Ward, Z. Wijs, Comp. Phys. Commun. 130, 260 (2000). 



11 



